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Laboratory experiments searching for galactic dark matter particles scattering off nuclei have 
so far not been able to establish a discovery. We use data from the XENONIOO experiment to 
search for dark matter interacting with electrons. With no evidence for a signal above the low 
background of our experiment, we exclude a variety of representative dark matter models that 
would induce electronic recoils. For axial-vector couplings to electrons, we exclude cross-sections 
above 6 x 10 -35 cm 2 for particle masses of m x = 2GeV/c 2 . Independent of the dark matter halo, 
we exclude leptophilic models as explanation for the long-standing DAMA/LIBRA signal, such as 
couplings to electrons through axial-vector interactions at a 4.4a confidence level, mirror dark matter 
at 3.6cr, and luminous dark matter at 4.6a. 

Dark Matter in the form of Weakly Interacting Mas¬ 
sive Particles (WIMPs) is typically expected to induce 
nuclear recoils in a terrestrial detector target [1] with 
an annually modulated rate due to the motion of the 
Earth around the Sun [2, 3]. Although such a modula¬ 
tion has been observed by the DAMA/LIBRA collabo¬ 
ration using sodium iodine [4], it is difficult to interpret 
it as a dark matter signal, given the null results from 
other experiments [5-7]. In fact, dark matter-induced 


nuclear recoils are excluded by these results unless one 
invokes models that are fine-tuned to create a signal only 
in DAMA/LIBRA but not in other experiments [8-10]. 
In contrast, dark matter-induced electronic recoils ap¬ 
pear as a viable explanation for the observed modula¬ 
tion because exclusions of other experiments do not ap¬ 
ply directly in this case [11-13]. We use data from the 
XENONIOO detector to rule out this possibility for three 
different, representative dark matter models. 
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FIG. 1: Fig. 1. Conceptual illustration of the anal¬ 
ysis. Shown is the DAMA/LIBRA rate (red) [20] with 
the modulated rate in (2 — 6) keV from the fit parameters 
in [4] (dark red). The distribution of the XENONIOO live 
time (blue) is indicated with its average background rate of 
5.3 events/(keV ■ tonne • day), which shows dents due to main¬ 
tenance or calibration campaigns. The region between the 
dashed lines (green) indicates the 70 summer live days where 
the modulated signal is expected to be largest. 


We interpret data from the XENONIOO detector that 
were acquired between February 28, 2011 and March 31, 
2012 for a total exposure of 224.6 live days and 34 kg fidu¬ 
cial mass. We have previously searched this data set for 
spin-independent [14] and spin-dependent [15] WIMP- 
induced nuclear recoils as well as for axion-induced elec¬ 
tronic recoils [16]. XENONIOO is located in the Gran 
Sasso underground laboratory. It consists of a liquid 
xenon target that is operated as a low-background time 
projection chamber [17]. Each particle interaction re¬ 
sults in two signals: The prompt scintillation signal 
(SI) is used here for energy estimation, and the de¬ 
layed ionization signal (S2) allows for 3D vertex recon¬ 
struction. Data reduction is performed in order to se¬ 
lect single-scatter low-energy (< lOkeV) recoils in the 
fiducial volume, while retaining maximal detector effi¬ 
ciency [16, 18]. At low energies, the remaining back¬ 
ground of XENONIOO is dominated by forward-scattered 
Compton events, resulting in a flat spectrum with a 
rate of 5.3 events/(keV • tonne • day) in the fiducial vol¬ 
ume [19] (File Al). This rate is more than two or¬ 
ders of magnitude lower than the average background 
rate of about 1019 events/(keV • tonne • day) reported by 
DAMA/LIBRA in the same energy interval [20, 21], and 
even smaller than their reported annual modulation am¬ 
plitude of (11.2 ± 1.2) events/(keV • tonne • day) [4]. Be¬ 
cause the DAMA/LIBRA collaboration has not pub¬ 
lished the composition of their background at low en¬ 
ergies, we test the minimum dark matter signal that 
would be required to cause the observed modulation. In 
this scenario, the constant spectrum is fully attributed 
to background, and only the modulated part itself is at¬ 
tributed to a 100% modulated dark matter signal as illus¬ 
trated in (Fig. 1). We ignore the practical difficulties of 


realizing such a highly modulated signal [3, 22] but con¬ 
servatively consider it as the case that is most challenging 
to exclude. The dark mat ter-induced rate would then be 
zero on December 2nd, and twice the measured modula¬ 
tion amplitude on June 2nd. It follows that there is an 
optimized time interval to consider for best sensitivity. 
To find this interval, the signal expected in XENONIOO 
was simulated for different time intervals centered around 
June 2nd. We take into account uncertainties from count¬ 
ing statistics in XENONIOO and DAMA/LIBRA, as well 
as the systematic uncertainty from the conversion of keV 
energy into SI [16]. The optimum time interval is found 
to be 70 live days around June 2nd, roughly correspond¬ 
ing to April 2011-August 2011 (Fig. 1) as indicated. Our 
expected sensitivity varies by less than O.lcr with changes 
of this interval of ±40 live days. A dedicated analysis of 
the time stability of XENONIOO electron recoil data will 
be presented elsewhere [23]. 

WIMP axial-vector coupling to electrons: A 
relativistic treatment of dark matter-electron scatter¬ 
ing shows that keV-scale electronic recoils can only be 
induced by dark matter particles with masses m x > 
lGeV/c 2 scattering inelastically off electrons with mo¬ 
menta on the order of MeV/c [11, 24]. A qualitatively sim¬ 
ilar result is obtained by a simple non-relativistic treat¬ 
ment of elastic two-body scattering. As shown in [11], 
even if the dark matter has tree-level (first-order) interac¬ 
tions only with leptons, loop-induced dark matter-hadron 
interactions dominate the experimental signatures and 
make the usual exclusions based on nuclear recoil analy¬ 
ses applicable. Thus, we consider here axial-vector A<g>A 
couplings between dark matter and leptons, since in this 
case, loop contributions vanish, while the WIMP-electron 
coupling is not suppressed by additional small factors of 
velocity v or mass ratio m e /m x . 
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FIG. 2: Fig. 2. Calculated ratio of the differential rates 

in xenon and sodium iodide for inelastic WIMP-electron scat¬ 
tering through axial-vector coupling. The structures around 
1 and 5 keV are owing to the small difference in the binding 
energies of the 3s and 2s shells in xenon and iodine. 


We use equation (30) in [11], with an additional fac¬ 
tor of 2 to account for electron occupancy from spin, to 
calculate the differential rate for WIMP-electron scatter¬ 
ing (File A2). The expected rate includes a sum over 
the atomic shells of the target, and for each shell, inte¬ 
grates the momentum wave function of the electrons to 
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get the contribution at a given recoil energy. Given the 
requirement that the energy deposited in the detector 
must be more than the binding energy of the electron, 
the largest contribution to the rate in a sodium iodide 
target comes from the 3s shell of iodine. The contribu¬ 
tions from sodium are two orders of magnitude smaller. 
The momentum-space wave functions for xenon atoms 
and iodine anions are nearly identical as a result of their 
similar electron structure. This has the important con¬ 
sequence that a comparison between sodium iodide and 
xenon is independent of the dark matter halo. The ratio 
of the calculated differential rates in xenon and sodium 
iodide are shown in Fig. 2 as a function of deposited en¬ 
ergy, considering the full shell structure. This ratio has 
negligible dependence on the WIMP mass. 

We contrast the DAMA/LIBRA signal, interpreted as 
WIMPs coupling to electrons through axial-vector inter¬ 
actions, with XENONIOO data. The energy spectrum of 
the modulation amplitude [4] is multiplied by the energy- 
dependent ratio from Fig. 2 and by a constant factor of 
1.88, which accounts for the time integral of the mod¬ 
ulated signal that is expected in our 70 summer live 
days (Fig. 1). The deposited electronic recoil energy in 
XENONIOO is estimated from the SI signal, measured 
in photoelectrons (PE), using the NESTv0.98 model [25] 
which consistently fits the available data [26-29]. The 
energy scale, shown in [16], includes a systematic un¬ 
certainty that decreases from 20% to 7% from 1 keV to 
10 keV, reflecting the spread and uncertainties in the 
measurements. The SI generation is modelled as a Pois¬ 
son process and the PMT resolution is taken into account 
in order to obtain the predicted XENONIOO SI spectrum 
from the scaled energy spectrum [18]. Our resolution is a 
factor 2 worse than that of DAMA/LIBRA; the feature 
at 5.2 keV in Fig. 2 is lost in this process. 

The converted DAMA/LIBRA and measured 
XENONIOO energy spectra are shown in Fig. 3. 
Part of the DAMA/LIBRA signal is expected to be 
seen below 2 keV due to the finite energy resolution of 
XENONIOO. The uncertainty in the converted signal 
includes both the statistical uncertainty in the original 
DAMA/LIBRA energy spectrum [4] as well as the 
uncertainties from our energy conversion. The electronic 
recoil cut acceptance, shown in [16], was applied to the 
converted DAMA/LIBRA spectrum. The uncertainty 
shown in the XENONIOO data is statistical. 

The energy region to determine the level of exclusion 
was chosen starting at the threshold of 3 PE [14] to 
the point where the DAMA/LIBRA signal falls below 
the expected average XENONIOO rate (cyan in Fig.3, 
calculated using a flat spectrum background model and 
scaled for the live time of the data set), which is at 
14 PE, corresponding to (2.0-5.9) keV. Taking system¬ 
atic uncertainties into account, a simple comparison of 
the integral counts in this energy interval excludes the 
DAMA/LIBRA signal as axial-vector coupling between 
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FIG. 3: Fig. 3. Contrasting XENONIOO data with 
DAMA/LIBRA. The DAMA/LIBRA modulated spectrum 
(red), interpreted as WIMPs scattering through axial-vector 
interactions, as it would be seen in the XENONIOO detector. 
The la band includes statistical and systematic uncertain¬ 
ties. The DAMA/LIBRA modulated spectrum interpreted as 
luminous dark matter is very similar, whereas the interpreta¬ 
tion as mirror dark matter is indicated separately (dark red). 
The (blue) data points are XENONIOO data from the 70 sum¬ 
mer live days with their statistical uncertainty. The expected 
average XENONIOO rate is also shown (dashed cyan). The 
shaded region from (3-14) PE was used to quantify the con¬ 
fidence level of exclusion. 


WIMPs and electrons at 4.4<r significance level, even con¬ 
sidering all events from the well-understood XENONIOO 
background [19] as signal candidates. To be consistent 
with previous analyses [16], the same data selection cuts 
were applied. The exclusion remains unchanged if we 
only impose a minimum set of requirements, namely that 
events have a single scatter in the fiducial volume with 
a prompt SI and delayed S2 signal in the correct energy 
range. Furthermore, the exclusion stays above 3 a confi¬ 
dence level even if we consider a 4.5cr downward deviation 
in the measured data points [26-28] that are used to set 
the energy scale, or if we set the light yield in xenon to 
zero below 2.9 keV, in contradiction with direct measure¬ 
ment [27, 28]. 

A profile likelihood analysis [30, 31] was performed to 
constrain the cross section cr 3 e = G 2 m 2 /n for WIMPs 
coupling to electrons through axial-vector interactions. 
To this end, we drop the assumption of a 100% mod¬ 
ulated rate and use the entire 224.6 live days data set. 
Fully analogous to [16], we use the same energy range and 
background likelihood function, derived from calibration 
data. We do not consider energy depositions below 1 keV, 
the lowest directly measured data point in [27]. The re¬ 
sulting XENONIOO exclusion limit (90% confidence level) 
is shown (Fig. 4) along with the lcr/2(7-sensitivity bands 
based on the background-only hypothesis. It excludes 
cross-sections above 6 x 10 -35 cm 2 for WIMPs with a 
mass of m x = 2GeV/c 2 . This is more than 5 orders of 
magnitude stronger than the one derived in [11] based on 
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FIG. 4: Fig. 4. Parameter space for WIMPs coupling 
to electrons through axial-vector interactions. The 

XENON100 upper limit (90% confidence level) is indicated 
by the blue line, along with the green/yellow bands indicat¬ 
ing the lcr/2cr sensitivity. For comparison, we also show the 
DAMA/LIBRA allowed region (red) and the constraint from 
Super-Kamiokande (SK) using neutrinos from the Sun, by 
assuming dark matter annihilation into rf or vv, both calcu¬ 
lated in [11]. 

data from the XENON10 detector, completely excludes 
the DAMA/LIBRA signal, and sets the strongest direct 
limit to date on the cross section of WIMPs coupling to 
electrons through axial-vector interactions.For compari¬ 
son, we also show the DAMA/LIBRA allowed region and 
the constraint from Super-Kamiokande using neutrinos 
from the Sun, by assuming dark matter annihilation into 
rf or vi>, both calculated in [11]. The XENON100 data 
completely excludes the DAMA/LIBRA signal and sets 
the strongest direct limit to date on the cross section of 
WIMPs coupling to electrons through axial-vector inter¬ 
actions, excluding cross-sections above 6 x 10 -35 cm 2 for 
WIMPs with a mass of m x = 2 GeV/c 2 . 

Kinematically Mixed Mirror Dark Matter: It 
has been suggested that multi-component models with 
light dark matter particles of ^MeV/c 2 mass might ex¬ 
plain the DAMA/LIBRA modulation [32]. A specific ex¬ 
ample of such a model, kinematically mixed mirror dark 
matter [33] , was shown to broadly have the right proper¬ 
ties to explain the DAMA/LIBRA signal via dark matter- 
electron scattering. In this model, dark matter halos are 
composed of a multi-component plasma of mirror parti¬ 
cles, each with the same mass as their standard model 
partners. The mirror sector is connected to the normal 
sector by kinetic mixing of photons and mirror photons at 
the level of ~ 10“ 9 , which provides a production mech¬ 
anism for mirror dark matter and a scattering channel 
with ordinary matter. While mirror hadrons would not 
induce nuclear recoils above threshold, mirror electrons 
(m' e = 511 keV/c 2 ) would have a velocity dispersion large 
enough to induce ~keV electronic recoils. 

The differential scattering rate of mirror electrons is 
proportional to gNn e i , where g is the number of loosely- 
bound electrons, assumed to be those with binding en¬ 


ergy < 1 keV [33], N is the number of target atoms 
and n e ’ is the mirror electron density.The detector- 
dependent quantities are N and g. In order to com¬ 
pare DAMA/LIBRA directly with XENON100, we ap¬ 
ply a constant scaling of gxe/ffNai • N Xe /N Nal = 0.89 to 
the DAMA/LIBRA spectrum and use the same proce¬ 
dure as in the case of axial-vector coupling: We again 
consider only the DAMA/LIBRA modulation signal, use 
the 70 summer live days, model scintillation in liquid 
xenon as described previously, and simply compare in¬ 
tegral counts up to the point where the DAMA/LIBRA 
signal falls below the expected average XENON100 back¬ 
ground data rate (at 13 PE), without background sub¬ 
traction. This excludes the DAMA/LIBRA signal as 
kinematically mixed mirror dark matter at 3.6cr confi¬ 
dence level. 

Luminous Dark Matter: The third model we con¬ 
sider is Luminous Dark Matter [34] , featuring a dark mat¬ 
ter particle with a ^keV mass splitting between states 
connected by a magnetic dipole moment operator. The 
dark matter particle upscatters in the Earth and later de- 
excites, possibly within a detector, with the emission of 
a real photon. The experimental signature of this model 
is a mono-energetic line from the de-excitation photon. 
A mass splitting S = 3.3 keV provides a good fit to the 
DAMA/LIBRA signal [34] which would be explained as 
scattering of a real photon from the de-excitation of a 
~GeV/c 2 dark matter particle that is heavy enough to 
undergo upscattering, but light enough to evade detec¬ 
tion in other direct searches. 

This signature is independent of the target material; 
only the sensitive volume affects the induced event rate. 
As rates are typically given per unit detector mass, scal¬ 
ing to volume is inversely proportional to target density. 
We thus apply a constant scaling factor to the differ¬ 
ential rate in DAMA/LIBRA which is the ratio of the 
target densities PNai/pxe = 1-29 in order to compare 
it to XENON100. Proceeding as in the previous two 
cases, we exclude the DAMA/LIBRA signal as luminous 
dark matter at 4.6u confidence level. Together with the 
other two exclusions presented above, this robustly rules 
out leptophilic dark matter interactions as cause for the 
DAMA/LIBRA signal. 
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